Non-monotonic pressure dependence of the dynamics of soft glass-formers at high 

compressions 
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In molecular dynamics simulations of soft glass-formers interacting via repulsions, we find that the 
glass transition temperature, dynamical heterogeneity, and glass fragility reach their maxima at the 
same crossover pressure Pd- Our analysis of the zero-temperature jammed states indicates that states 
at Pd have the highest bond orientational order with the largest spatial fluctuation. Correspondingly, 
the low-frequency normal modes of vibration are the least localized and the average potential energy 
barrier along these modes are the highest for jammed states in the vicinity of Pd- The reentrant 
glass transition and dynamics of supercooled liquids are thus predictable by these structural and 
vibrational precursors in the zero-temperature jammed states. 
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I. INTRODUCTION 

Upon fast cooling or compression, a liquid or colloidal 
suspension undergoes the glass transition with sluggish 
dynamics. Due to the lack of convincing phase transi- 
tion signatures, the nature of the glass transition is still 
controversial [lPrlG| . In practice, the glass transition tem- 
perature is usually defined as the temperature at which 
the relaxation time exceeds the measurable time window 
or is extrapolated from the divergence of a functional 
fit, e.g. the Vogel-Fulcher-Tamman (VFT) or mode 
coupling [BJ form, to the relaxation time data. Previous 
studies have shown that the glass transition tempera- 
ture usually increases monotonically with the pressure 
or equivalent ly the volume fraction, so does the glass 
fragility |13 - [l4| . Only very recently, it has been noticed 
that such a picture may be violated by core-softened col- 
loids at high pressures. For instance, a recent simulation 
has shown that the mode coupling glass transition tem- 
perature of a highly compressed model glass-former is 
reentrant as a function of the volume fraction [l5[ . 

Core-softened colloids can exhibit complicated and rich 
phase behaviors and dynamics at high pressures [H|-|lQj], 
which have been mostly attributed to the softness of the 
colloids. However, the underlying mechanism is appar- 
ently lacking. Furthermore, our understanding of the dy- 
namics of soft glass-formers at high compressions is still 
poor. To our knowledge, the reentrant glass transition 
temperature [l5| is the only phenomenon ever reported 
in the literature. Other interesting and important as- 
pects of soft glass-formers at high pressures, e.g. dynam- 
ical heterogeneity and glass fragility, have not yet been 
concerned about. Approaching the glass transition, the 
dynamics of sup ercooled liquids are spatially heteroge- 
neous 0, EH [Hrllll • The dynamical heterogeneity tends 
to grow near the glass transition, implying a diverging 
dynamic correlation length. As one of the most striking 
dynamical anomalies, will the dynamical heterogeneity 
show non-monotonic pressure dependence as well? 

Near the glass transition a supercooled liquid loses er- 



godicity, caged in local potential energy basins over a 
long time before being able to escape. From this point 
of view, the slow dynamics and dynamical heterogeneity 
of supercooled liquids are determined by the local fluc- 
tuation of the structural softness and the quasi-localized 
nature of the low-frequency normal modes of vibration 
of the metastableglass states at the local potential en- 
ergy minima [26|-|30(. Due to the structural disorder, the 
low frequency modes of the metastable glasses are quasi- 
localized, i.e. all particles are involved in the vibration 
but a small fraction in localized regions vibrate more 
strongly [30| - |33| . The quasi- localization is a special fea- 
ture of the low-frequency vibrations of disordered solids 
and different from the Anderson localization in which 
most of the particles stay still except for a localized re- 
gion. The low-frequency vibrations are usually localized 
in soft spots, i.e. more disordered regions of loosely con- 
strained particles }27l - l30j . 

A recent study of the zcro-tcmperature (T = 0) 
marginally jammed solids has shown that along more 
quasi-localized low-frequency modes the potential energy 
barrier heights are lower [32[. Here jamming is restricted 
to packings of frictionless spheres interacting via repul- 
sions at T = 0. A packing is jammed when it be- 
comes rigid with nonzero elastic moduli and a coordi- 
nation number larger than the isostatic value, i.e. the 
minimum number of constraints per particle to maintain 
global mechanical stability j34|. The jammed states are 
actually metastable glasses with purely repulsive interac- 
tions. The observed correlation between the energy bar- 
rier height and quasi- localization of jammed states sug- 
gests that a metastable glass with more quasi-localized 
low-frequency vibrations would be more vulnerable to ex- 
citations and consequently have a lower glass transition 
temperature. To our knowledge, however, this picture 
has not yet been quantitatively demonstrated. 

The study reported in this paper is inspired by 
our recent observation that at T = there exists 
a crossover volume fraction (pd (or correspondingly a 
crossover pressure Pd) that separates marginally jammed 
solids from deeply jammed ones with distinct proper- 
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ties [19|. Our analysis indicates that the low- frequency 
quasi-localization gets weaker with increasing the vol- 
ume fraction for marginally jammed solids until (f>d is 
attained, after which the quasi-localization of deeply 
jammed solids exhibits opposite volume fraction depen- 
dence [Tj|. From the picture discussed above, the reen- 
trant quasi-localization at 4>d must have nontrivial effects 
on the glass transition and dynamics of supercooled liq- 
uids. 

In this paper, we study the dynamics and structure of 
supercooled liquids consisting of soft particles in a wide 
range of pressures via molecular dynamics (MD) simu- 
lations. We find that not only the glass transition tem- 
perature but also the dynamical heterogeneity and glass 
fragility are all reentrant near the crossover pressure Pd- 
We propose an explanation of this non-monotonic pres- 
sure dependence from the analysis of the corresponding 
T = jammed states. Our study provides quantitative 
evidence of the direct link between the glass transition 
temperature and low-frequency quasi-localization of the 
normal modes of the T = metastable states. 



II. SIMULATION DETAILS 

Our model systems are three dimensional cubic boxes 
with periodic boundary conditions filled with a binary 
mixture of N = 1000 frictionless spheres with the same 
mass m. We also study some N = 10000 systems and 
do not observe significant finite size effects. The diame- 
ter ratio of the large to small spheres is 1.4 to effectively 
avoid crystallization. The potential between two parti- 
cles i and j is V (rij) = e (1 — m / <Jij ) /a when their 
separation is smaller than the sum of their radii <xy , 
and zero otherwise, a is a tunable parameter to deter- 
mine the softness of the interaction. The same model 
system has been widely applied to the study of jamming 
and glass dynamics. We study both harmonic (a = 2) 
and Hertzian (a = 5/2) repulsions. In this paper, we 
only present the results for harmonic repulsion. Our ma- 
jor findings are valid for Hertzian repulsion as well. We 
set the characteristic energy scale e, small particle diam- 
eter cr, and particle mass m to be the units. The time 
and temperature are in the units of a/y^e/m and e/ks, 
respectively, where ks is the Boltzmann constant. 

For glass-formers at T > 0, we perform MD simula- 
tions at constant temperature and pressure to obtain the 
time evolution of particle motion according to the equa- 
tions of motion 1351: 



dr; _ 

— = vi + \n, 

at 

dt m 



n ^ — ' 



(i) 

(2) 
(3) 



Fij = — VVij is the force acting on particle i by particle j, 
L is the length of the simulation box, and £ and A are La- 
grange multipliers to maintain constant temperature and 
pressure. We apply 4-variable Gear predictor-corrector 
algorithm to integrate the equations numerically (35| . 

The relaxation time of the liquids is measured from the 
self-part of the intermediate scattering function 

f s (k, t) = | ex p(^ • (*) - ^ w 



where the sum is over all large particles, rj(i) is the lo- 
cation of particle j at time t, and k is chosen in the 
a;— direction with k = \k\ satisfying the periodic bound- 
ary conditions and being approximately the value at the 
first peak of the static structure factor. The relaxation 



time r is determined by F s (k,r) 



l F s (k,0). All the 



measures at T > discussed in this paper are taken after 
the system has been equilibrated for several r. 

Widely accepted tools to probe the dynamical hetero- 
geneity of supercooled liquids include the non-Gaussian 
parameter a 2 [2l| . _36| ._ 37| and four-point dynamical sus- 



ceptibility X4 
defined as 



The non-Gaussian parameter is 



a 2 {t) 



(Ar(t) 4 ) 



(1 + 2/d) (Ar(t) 2 Y 



(5) 



where Ar(t) = r(t) — r(0) is the particle displacement 
at time t, d is the dimension of space, and (.) denotes 
the time average over all the particles. For simple liq- 
uids satisfying the Gaussian distribution, a 2 = 0. For 
supercooled liquids, however, a 2 is greater than zero and 
shows interesting time and temperature dependence [2l[ . 
The four-point dynamical susceptibility %4 captures the 
fluctuation of the number of mobile particles [3J, l22|, |23J . 
It measures the variation of an overlap function Q(a,t): 



X4(a,t) = -[(Q(a,t) 2 ) - (Q(a 7 t)Y 



(6) 



where (.} denotes the time average. The overlap function 
Q(a,t) = J2iLi W a (\ri(t) — f^(0) |) evaluates the similar- 
ity between two configuration snapshots separated by a 
time interval t, where a is a preset length, and W a (x) = 1 
if x < a and zero otherwise. Note that here we only con- 
sider the self part of the overlap function, which domi- 
nates the inter-particle parts and describes the dynamical 
heterogeneity well [13, |23| . Both a 2 and %4 are larger if 
the dynamics arc more heterogeneous. 

The structure of the glass-formers is evaluated from 
the pair distribution function of large particles j35| 



g(r) 



{N/2f 



EE* 



(7) 



where and vi arc the position and velocity of particle i, 



where the sums are over all the large particles and {.) 
denotes the time average. 
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We obtain jammed configurations at T — by quickly 
quenching ideal gas states to local potential energy min- 
ima using L-BFGS method [Hj. We tune the volume 
fraction or equivalently the particle size successively un- 
til a desired pressure is attained. At each pressure, we 
generate over 1000 distinct jammed states and take the 
average over them. 

In the T = jammed states, the local bond orienta- 
tional order of particle i is defined as [H, 



>(0 = 



\ 



(8) 



where 



Q&m(i) 



m— — 6 



N b (i) 



3 = 1 



(9) 



with Nb(i) the number of nearest neighbors of particle i 
determined by the Voronoi tessellation, and Y§ m (fij ) the 
spherical harmonics. We measure the average bond ori- 
cntational order (Qe) and its spatial fluctuation SQq = 



(Qe) ~ (Qe) j where (.) denotes the particle and con- 
figuration average. 

We diagonalize the Hessian matrix of the T = 
jammed states using ARPACK [4l[ to obtain the nor- 
mal modes of vibration. For each mode, we calculate its 
participation ratio 
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FIG. 1: Self-part of the intermediate scattering function 
F a (k,t) with || = 10 measured at T = 0.003 and various 
pressures for systems with harmonic repulsion (a = 2). The 
pressures are 0.04 (black solid), 0.07 (red dotted), 0.08 (blue 
dashed), and 0.09 (green dot-dashed) in (a) and 0.2 (black 
solid), 0.26 (red dotted), 0.3 (blue dashed), and 0.6 (green 
dot-dashed) in (b). 



with Rq and the energy minimization will lead state R 
back to i?o- We thus define the potential energy bar- 
rier height along mode n, AVJnax(w n ) as the maximum 
AV above which state R is no longer in the same basin 
with Rq. Apparently, more stable jammed states should 
have higher AV max and hence possibly higher glass tran- 
sition temperatures. In the following sections we will 
show the interesting correlations between the measure of 
Qe, AV^ lax , p, £, and the glass transition temperature. 



Ei=l \ e r 



(10) 



where uj n and e n ^ are the frequency of the n th mode and 
polarization vector of particle i in the mode. The par- 
ticipation ratio measures the extensiveness of a mode, 
i.e. the fraction of particles effectively involved in the vi- 
bration. More localized modes have smaller participation 
ratios and vice versa. We also measure the correlation 
function of the polarization vectors in mode n [l|| : 

r m - ^ J=l e "'' ' e "' J ^ r ~ n ^ fin 

n{ ) ~ T N T N S(r-r- ) ' ( j 

The normalized correlation function C„ (r) = 
C n (r)/C n (Q) can be well fitted with 



C^(r)=C n0 exp(-r/0+A J 



(12) 



where C„o and A are fitting parameters [19f. A more 
localized mode has a smaller correlation length £. 

We use i?o to denote a T = jammed state in 
the configurational space. If we move the state along 
mode n to R = Rq + ue n , the potential energy rises to 
V{R) = V(Ro) + AV. If |u| is small enough, the per- 
turbed state R is still in the same potential energy basin 



III. DYNAMICS AND STRUCTURE OF 
SUPERCOOLED SOFT GLASS-FORMERS 



Fig. Q] shows the pressure evolution of the self-part of 
the intermediate scattering function F s (k,t) at a fixed 
temperature T = 0.003 for systems with harmonic re- 
pulsion, which simulates the route of the colloidal glass 
transition. As shown in Fig. [TJa), the relaxation time 
increases rapidly with increasing the pressure as usually 
observed d, [H, Ell- When the pressure gets higher, a 
typical two-step relaxation emerges due to the cage effect 
and becomes more and more pronounced. From our pre- 
vious experience, we may naturally expect that the glass 
transition happens at a critical pressure above which the 
systems turn into glasses, which is true for hard spheres 
[42| or soft spheres at low enough temperatures 0, [l]| . 
At T = 0.003, however, we surprisingly find that above a 
crossover pressure the relaxation time decreases with fur- 
ther compression. As shown in Fig. [ljb), the dynamics 
get faster at high compressions and the two-step relax- 
ation eventually disappears at high enough pressures. 

In order to have a quantitative picture of the pressure 
dependence of the glass transition, we measure the re- 
laxation time in equilibrium at various temperatures and 
pressures. At constant pressure, we fit the relaxation 







FIG. 2: (a) — (b) Arrhenius plots of the relaxation time t, 
(c) VFT glass transition temperature To at which r diverges 
(black circles), glass transition temperature T g at which r = 
10 6 (red squares), and (d) glass fragility k for systems with 
harmonic repulsion (a = 2). The pressures are 0.00001 (black 
circles), 0.00005 (red squares), 0.0005 (blue diamonds), 0.005 
(green triangles), 0.05 (orange stars), and 0.1 (violet pluses) 
in (a), and 0.2 (black circles), 0.25 (red squares), 0.3 (blue 
diamonds), and 0.4 (green triangles) in (b). The solid curves 
are the fits with Eq. (Tl3)) . 

time with the VFT function 

T = T CXp ( — , (13) 

where To, A, and Tq are fitting parameters. To is the 
VFT glass transition temperature at which r = oo. 

In panels (a) and (b) of Fig. [5] we plot the relaxation 
time measured at constant pressure as a function of T g /T , 
i.e. the Arrhenius plot, where T g is defined as the tem- 
perature at which r = 10 6 . At all the pressures, our data 
can be well fitted with Eq. ([13]) , from which we can esti- 
mate To and T g . As shown in Fig. &[c), both Tq and T g 
vary non-monotonically with the pressure, in agreement 
with the recent observation that the mode coupling glass 
transition temperature is reentrant as a function of the 
volume fraction [l5j . Tq and T g reach their maximum 
values (~ 0.003) at a crossover pressure Pd w 0.2, which 
is approximately equal to the critical pressure separat- 
ing marginal jamming from deep jamming at T = fl9j |. 
This agreement is not just an coincidence, which will be 
discussed in detail in Section ITVl 

The slope of the Arrhenius plot at T = T g , k = 

^j^Jt) It=t 8 measures the glass fragility. A fragile 
(strong) glass has a large (small) k. Fig. f^a) shows that 
the glass fragility increases with the compression when 
the pressure P < Pd, in consistent with previous obser- 
vations [TH, [l3j]. However, this trend holds only up to 
P ps Pd, as shown in Fig. H^b) that the glasses become 
stronger instead when P > Pd- The explicit pressure de- 
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FIG. 3: Dynamical heterogeneity of supercooled liquids with 
harmonic repulsion (a = 2) measured at various pressures and 
a fixed temperature T = 0.003. Panels (a) and (b) show the 
non-Gaussian parameter ot2{t). The pressures are 0.04 (black 
solid), 0.07 (red dotted), 0.08 (blue dashed), and 0.09 (green 
dot-dashed) in (a) and 0.2 (black solid), 0.26 (red dotted), 
0.3 (blue dashed), and 0.6 (green dot-dashed) in (b). Panels 
(c) — (f) are the contour plots of the four-point dynamical sus- 
ceptibility xt( a ,t) measured at pressures of (c) 0.06, (d) 0.08, 
(e) 0.4, and (d) 0.6. 



pendence of k is plotted in Fig. EJd) with a maximum 
present at Pd, indicating that glasses at Pd are the most 
fragile. The gap between Tq and T g also reflects the glass 
fragility to some extent, which can be seen from the ex- 
pression of k assuming that the relaxation time satis- 
fies Eq. (fHij) . Strong glasses may have large values of 
T g — Tq due to their more Arrhenius behaviors, which is 
exactly the case in Fig. Htc). It is interesting to know 
that high pressure has nontrivial effects on glass proper- 
ties: it can melt glasses and make glasses stronger. From 
panels (c) and (d) of Fig. [2] we can see that with com- 
parable fragility glasses at P > Pd can have much higher 
glass transition temperature than those at P < Pd- Our 
observation here thus proposes an effective way to obtain 
strong glasses with high glass transition temperatures. 

Accompanied with the drastic slowdown upon the glass 
transition, the dynamics of supercooled liquids are spa- 
tially heterogeneous. The dynamical heteroge neity grows 
with the increase of the relaxation time [H, [22|,[23[. Since 
the relaxation time is reentrant in pressure at constant 
temperature, is it possible that the dynamical hetero- 
geneity exhibits the similar pressure dependence? The 
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FIG. 4: (a) Pair distribution function of large particles g(r) 
of supercooled liquids with harmonic repulsion (a = 2) mea- 
sured at a fixed temperature T = 0.003 and pressures of 0.01 
(black solid), 0.06 (red dotted), 0.2 (blue dashed), 0.4 (green 
dot-dashed), and 0.5 (violet dot-dot-dashed), (b) Crossover 
pressure Pj at which the first peak of g(r) reaches the maxi- 
mum height at constant temperature. The red solid line shows 
the scaling relation Pj/4> 2 ~ T 1//2 which fits the data well up 
to the maximum glass transition temperature T™ ax ps 0.003. 



answer is positive from Fig. [3] which shows the non- 
Gaussian parameter a 2 and four-point dynamical suscep- 
tibility \4 measured at various pressures on both sides of 
Pd and a fixed temperature T = 0.003. 

At all the pressures the non-Gaussian parameter a 2 (t) 
is peaked at a time t a scaled with the relaxation time 
t. The magnitude of the peak signifies the strength of 
the dynamical heterogeneity. Panels (a) and (b) of Fig. [3] 
indicate that a 2 behaves opposite pressure dependence 
on the two sides of Pd- When P < Pd, the peak of a 2 
grows at an increscent t a , meaning the growth of the 
dynamical heterogeneity upon compression. When P > 
Pd, however, the peak of a 2 drops instead and moves to 
a shorter t a . 

Panels (c)-(f) of Fig. [3J are contour plots of the four- 
point dynamical susceptibility Xi{ a it) measured at dif- 
ferent pressures. \i reaches its maximum at a length 
a « 0.4 weakly depending on the pressure and a time t x 
showing similar pressure dependence as r and t a . Along 
with the compression, the peak value of X4 initially grows 
until P w Pd and then drops, as what a 2 does. Like a 2 , 
large %4 means strong dynamical heterogeneity. Fig. [3] 
thus indicates that at the crossover pressure Pd, together 
with the emergence of the maximum glass transition tem- 
perature and glass fragility, the dynamics are the most 
heterogeneous. 

Our recent study of the T = jammed states suggests 
that the crossover pressure Pd may have its structural 
origins [19} . At Pd, the coordination number of jammed 
states is 12, implying that particles may start to inter- 
act with their second nearest neighbors. Structurally a 
second peak thus emerges at r < ctl in the pair distri- 
bution function of large particles g(r), where <tl is the 
diameter of large particles. For soft glass formers stud- 
ied here, we observe the similar structural change across 
Pd- Fig. HJa) shows the pressure evolution of g(r) at 



T = 0.003. Compared to the T = results the ther- 
mal motion smears out the discontinuous jump of g(r) 
at r = a i, and masks the initial formation of the second 
peak on the left hand side of r = <tl at P w Pd- However, 
the emergence of the second peak (or even more peaks 
at extremely high pressures) at r < ctl is still a robust 
signature that distinguishes the glass-formers on the two 
sides of Pd- 

It has been shown that during the formation of soft 
colloidal glasses at a fixed low temperature the first peak 
of g(r) reaches the maximum height at a crossover pres- 
sure Pj , which is embedded in the glass regime and remi- 
nisces the T = jamming transition 0, 0, |44|. Beyond 
the maximum glass transition temperature T™ ax ~ 0.003 
(sec Fig. because there are no solid states, Pj is no 
longer associated with jamming. As shown in Fig. Hljb), 
at low temperatures up to T™ ax , Pj/<fi 2 is scaled well 
with T(«-!)/« as reported HI Si- When T > T^ ax , 
it is interesting that the scaling breaks down and Pj re- 
mains almost constant in the temperature, corresponding 
to the decrease of the volume fraction [45[ . 

All the results discussed in this section indicate that 
the existence of the crossover at Pd is not an accident, 
because all the quantities that we concern about undergo 
apparent changes there. According to the proposal that 
the dynamics of supercooled liquids reflects the structural 
and vibrational properties of the T = metastable states 
[26l - [30l ] , we may be able to seek the origin of the reentrant 
glass transition from the analysis of the T = jammed 
states. 



IV. IMPLICATIONS FROM T = JAMMED 
STATES 

We recently find that there exists a crossover volume 
fraction cf>d « 1.2 that divides the T = jammed states 
into marginally and deeply jammed regimes with distinct 
structures and power-law scaling of typical quantities 
such as the potential energy, elastic moduli, and coordi- 
nation number [l9j . The pressure at the crossover agrees 
with Pd- We have also made a prediction that the non- 
monotonic pressure dependence of the quasi-localization 
of the low-frequency normal modes of vibration leads to 
the reentrant glass transition. In this section, we will 
show that the unusual dynamics at high compressions 
discussed in section [TTT] are indeed strongly correlated to 
the properties of the T = jammed states. 

Upon compression, each particle interacts with more 
and more neighbors. It thus sounds plausible that the 
heterogeneity induced by the particle size dispersion may 
be suppressed and the jammed states may be more 
and more ordered in structure. Panels (a) and (b) of 
Fig. [5] show the potential energy distribution in a two- 
dimensional cross section of the jammed states on both 
sides of Pd, scanned by a point tracer. When P < Pd, we 
can clearly see the boundary of each particle, indicating 
that particles only interact with their nearest neighbors. 
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FIG. 5: Structure of the T — jammed states with harmonic 
repulsion (a = 2). Panels (a) and (b) are the potential energy 
field in a two-dimensional cross section of the states at P = 
0.00001 and 0.5, respectively. Panels (c) and (d) show the 
pressure dependence of the average bond orientational order 
(Qe) and its relative fluctuation SQs/ (Qe). 




When P > Pd, it is hard to identify single particles due 
to the penetration of the particle interactions to farther 
neighbors. The potential energy field is still heteroge- 
neous in space with no visible increase of the structural 
order. 

We then show in Fig. [5jc) the bond-orientational order 
(Qe) of jammed states as a function of the pressure. A 
state with a larger (Qe) is more ordered and vice versa. 
When P < Pd the jammed states are more ordered upon 
compression, as expected. It is interesting that the order 
stops increasing at Pd- When P > Pd, the compres- 
sion boosts the disorder instead. Recent studies have 
demonstrated that the particle displacement in a super- 
cooled liquid is stro ngly correlated to the structural order 
of the particle [27], l30| . Disordered particles with small 
Qe are usually loosely constrained and tend to move to 
longer distances than ordered particles. Therefore, dis- 
ordered states should be more vulnerable than ordered 
states to the same strength of external excitations. From 
this phcnomcnological picture, it is reasonable to expect 
a reentrant glass transition around Pd, which is exactly 
observed in our MD simulations. 

The fluctuation of the bond-oricntational order, 
5Qe/ (Qe) behaves similarly to (Qe), as shown in 
Fig. EJd). The fluctuation of the local structural or- 
der reaches the maximum at Pd. In consideration of 
the correlation between the particle dynamics of super- 
cooled liquids and structural order of metastable states, 
the reentrant dynamical heterogeneity in the vicinity of 
Pd originates from the pressure dependence of the struc- 
tural order fluctuation of the T = jammed states. 

The normal modes of vibration are the fundamentals 
to understand the properties of solids under excitations, 
e.g. the energy transport. Recent studies have suggested 



FIG. 6: Pressure dependence of the participation ratio (p), 
correlation length of the polarization vectors {£}, and poten- 
tial energy barrier height along vibrational modes (AV ma x) 
averaged over the 20 lowest frequency normal modes of vi- 
bration for T — jammed states with harmonic repulsion 
(a — 2). The inset to (b) shows the normalized correlation 
function C%(r) at P = 0.001 (black circles) and 0.2 (red 
squares). The lines are the fits with Eq. (|12|l from which £ is 
extracted. 



that the low-frequency modes play important roles in the 
determination of the heterogeneous dynamics of super- 
cooled liquids or sheared glasses |26l - [3fj |. As discussed 
above, disordered particles with small Qe are loosely con- 
strained and easy to move. In the low-frequency vibra- 
tions, these particles also form local soft spots with longer 
polarization vectors than the others and result in the 
quasi-localization. We may then be able to find precur- 
sors of the reentrant glass transition in the analysis of 
the low-frequency modes. 

Fig. |H|a) shows the average participation ratio (p) of 
the 20 lowest frequency modes, (p) reaches the maximum 
at Pd, meaning that the low- frequency modes at Pd are 
the least localized, which is consistent with the presence 
of the maximum order. To quantitatively illustrate the 
pressure dependence of the localization, we extract the 
correlation length £ from the correlation function of the 
polarization vectors C n (r) (see Section [IT] for details). In- 
terestingly, Fig. [^b) shows that £ is the largest at Pd as 
well, which is another strong evidence of the least local- 
ization at Pd [l9j . 

In the potential energy landscape, each jammed state 
sits at the local potential energy minimum of a basin 
of attraction. When thermally excited, the system ex- 
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FIG. 7: Participation ratio p(w) of two three-dimensional T = 
metastable states consisting of 1000 particles interacting 
via Lennard- Jones (black circles) and repulsive Lennard- Jones 
(red squares) interactions at a reduced number density p = 
1.2. The systems are exactly the same as in Ref. 46. 

plores the nearby configurational space. Each basin has 
a potential energy barrier in any direction. If the kinetic 
energy is large enough, the system will overcome the bar- 
riers and the glass melts. In this picture, the glass tran- 
sition temperature is correlated to the potential energy 
barrier heights. In Fig. [6jc) we show the potential en- 
ergy barrier height {Ay max ) averaged over the 20 lowest 
frequency modes. The jammed states around Pd are the 
most stable with approximately the highest (AVi nax ), al- 
though the peak of (AVm ax ) is not exactly at Pd probably 
due to some uncertainties in the calculation such as the 
large configurational variation of AVmax and unphysical 
effects of the energy minimization algorithm. 

From the analysis of the structural order and normal 
modes of vibration of the T = jammed states, we are 
able to understand the non-monotonic pressure depen- 
dence of the dynamics at high compressions. Here we 
show the strong couplings between the structural order, 
quasi- localization, and energy barrier height. All these 
quantities are reentrant around Pd, acting as the precur- 
sors of the reentrant glass transition. 

V. CONCLUSIONS 

In this study, we find that soft glass-formers with 
purely repulsive interactions behave non-monotonic dy- 
namics upon compression. There exists a crossover pres- 
sure Pd at which the glass transition temperature, glass 
fragility, and dynamical heterogeneity reach their max- 
ima. Although a large number of studies have been per- 
formed to understand the elusive glass transition prob- 
lems, the reentrant glass transition and dynamics are 
still quite unexpected and may thus raise new challeng- 
ing questions. The unusual pressure dependence of the 
highly compressed glasses may have some potential ap- 



plications. For instance, we can produce strong glasses 
with high glass transition temperatures using highly com- 
pressed core-softened colloids. Although strong glasses 
can be obtained at low compressions as well, the glass 
transition temperature has to be extremely low as the 
compensation (see Fig. [2|. 

The results reported in this paper require experimental 
verification. Possible candidates include extremely soft 
colloids such as star polymers, charged colloids with long 
range electrostatic interactions, or dusty plasmas. Here 
we only concern about purely repulsive systems. It is 
also interesting to know whether and how the attraction 
alters the picture. 

The crossover pressure Pd is consistent with that sep- 
arating marginal jamming from deep jamming at T = 
p^ |. which inspires us to search for origins of the reen- 
trant glass transition from the analysis of the T = 
jammed states. Interestingly, we observe strong correla- 
tions between the structural order and properties of the 
normal modes of vibration of the T = jammed states. 
The jammed states at Pd are the most stable with the 
highest bond oricntational order, the largest spatial fluc- 
tuation of the order, the least localized low-frequency 
modes, and the highest potential energy barriers. It is 
thus not hard to understand why the glass transition is 
the highest and the dynamics are the most heterogeneous 
at P d . 

Due to the strong correlation between the glass 
transition temperature and properties of metastable 
glasses such as the structural order, low-frequency quasi- 
localization, and energy barrier height, we are able to 
predict and understand the effects of some conditional 
changes (e.g. type of interactions and boundary condi- 
tions) on the glass transition from the analysis of the 
metastable states, e.g. the participation ratio. For in- 
stance, recent studies have shown that the inclusion of 
the attraction significantly influences the dynamics of su- 
percooled liquids: A Lennard- Jones (LJ) system has a 
higher glass transition temperature than its counterpart 
with repulsive Lennard- Jones (RLJ) interactions (46ll47||. 
This result is predictable from the participation ratio of 
the low-frequency vibrational modes of the metastable 
states. The low-frequency modes of the RLJ system must 
be more localized than those of the LJ systems, which is 
exactly the case shown in Fig. [7] We believe that the 
participation ratio of low-frequency vibrational modes of 
metastable states is a simple but efficient tool to predict 
the behaviors of supercooled liquids in various systems. 
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